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1.  Summary  of  research 

We  have  performed  research  on  the  algorithm  development,  analysis,  implemen¬ 
tation  and  application  for  high  order  weighted  essentially  non-oscillatory  (WENO) 
finite  difference  and  finite  volume  schemes  and  related  methods,  for  solving  computa¬ 
tional  fluid  dynamics  (CFD)  problems  and  other  applications  containing  strong  shock 
waves  and  complicated  smooth  region  structures.  The  goal  of  this  research  is  to  help 
obtaining  more  robust,  cost  effective,  and  reliable  numerical  tools  for  computational 
fluid  dynamics  problems  and  other  applications,  that  are  of  interest  to  AFOSR. 

There  are  22  publications  in  refereed  journals,  which  have  appeared  or  have  been 
accepted  for  publications  and  have  quoted  partial  support  by  this  grant,  see  section 

2. 

In  [11]  (the  numbers  here  and  below  refer  to  that  of  publications  listed  in  sec¬ 
tion  2),  we  have  developed  a  fifth  order  weighted  essentially  non-oscillatory  (WENO) 
fast  sweeping  method  for  solving  static  Hamilton-Jacobi  equations,  utilizing  two  ap¬ 
proaches  of  accurate  inflow  boundary  condition  treatment  which  allows  the  usage  of 
Cartesian  meshes  regardless  of  the  domain  boundary  shape.  This  method  is  expected 
to  very  useful  in  a  wide  range  of  applications  involving  level  set  methods  and  front 
propagation. 

In  [1],  we  have  developed  high  order  essentially  non-oscillatory  (ENO)  Lagrangian 
schemes  based  on  the  Lax-Wendroff  (LW)  type  time  discretization  procedure.  Exten¬ 
sive  numerical  examples  are  presented,  for  both  the  second  order  spatial  discretiza¬ 
tion  using  quadrilateral  meshes  and  third  order  spatial  discretization  using  curvilinear 
meshes.  These  methods  are  expected  to  be  very  useful  for  multi-phase  flows  and  they 
can  capture  phase  boundaries  sharply. 

In  [15],  we  have  developed  a  framework  to  construct  uniformly  high  order  accu¬ 
rate,  total  variation  diminishing  finite  volume  and  discontinuous  Galerkin  schemes 
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for  solving  one  dimensional  nonlinear  scalar  conservation  laws.  The  main  idea  is  to 
measure  the  total  variation  of  the  complete  piecewise  polynomial,  to  design  a  simple 
scalar  limiter  which  keeps  high  order  accuracy,  and  to  follow  exact  time  evolution. 
The  method  in  this  paper  overcomes  the  traditional  barriers  of  second  order  accuracy 
for  total  variation  diminishing  schemes.  However,  the  requirement  of  exact  time  evo¬ 
lution  prevents  its  generalization  to  multi-dimensions  and  systems.  In  [16],  we  have 
adapted  and  generalized  this  framework  to  construct  uniformly  high  order  accurate  fi¬ 
nite  volume  and  discontinuous  Galerkin  schemes  satisfying  a  strict  maximum  principle 
for  scalar  conservation  laws.  The  main  breakthrough  includes  a  framework  allowing 
the  maximum  principle  for  the  cell  averages  with  Euler  forward  time  discretization 
without  affecting  spatial  high  order  accuracy,  and  a  simple  scaling  limiter  which 
requires  the  evaluation  of  the  polynomial  solution  only  at  fixed  quadrature  points 
and  keeps  high  order  accuracy.  This  breakthrough  allows  us  to  design  high  order 
schemes  satisfying  a  strict  maximum  principle  for  multi- dimensional  scalar  problems 
and  passive  convection  in  a  divergence- free  velocity  field,  including  two-dimensional 
incompressible  Euler  equations  in  vorticity-streamfunction  formulation.  The  method¬ 
ology  is  extended  in  [17]  to  Euler  systems  of  compressible  gas  dynamics  to  maintain 
positivity  for  density  and  pressure  without  affecting  high  order  accuracy,  and  in  [18] 
to  compressible  Euler  equations  with  geometric,  reaction,  gravity  and  other  types 
of  source  terms.  An  adaptation  of  the  methodology  to  high  order  finite  difference 
WENO  schemes  is  carried  out  in  [20],  and  a  minimum  entropy  principle  of  high  order 
schemes  for  gas  dynamics  equations  in  [21].  Finally,  generalizations  to  arbitrary  un¬ 
structured  triangular  meshes  are  carried  out  in  [22],  This  line  of  research  has  led  to 
the  first  genuinely  uniformly  high  order  schemes  which  satisfy  strict  maximum  prin¬ 
ciple  and/or  positivity-preserving  properties  as  appropriate,  which  are  very  simple  to 
implement  and  have  shown  excellent  performance  in  very  demanding  numerical  tests 
such  as  detonation  waves  in  [8].  We  expect  this  methodology  to  be  widely  used  in 
applications.  An  extension  to  convection  diffusion  equations  for  a  non-standard  finite 
volume  scheme  is  established  in  [14].  The  student  co-author,  Dr.  Xiangxiong  Zhang 
(who  obtained  his  Ph.D.  degree  in  2011),  has  been  chosen  as  one  of  the  three  winners 
in  the  prestigious  SIAM  Student  Paper  Competition  for  his  joint  paper  [16]. 

In  [5],  we  have  developed  a  high  order  finite  difference  numerical  boundary  con¬ 
dition  for  solving  hyperbolic  conservation  laws  on  a  Cartesian  mesh.  The  challenge 
results  from  the  wide  stencil  of  the  interior  high  order  scheme  and  the  fact  that  the 
physical  boundary  may  not  be  aligned  with  the  grid  points  in  the  Cartesian  mesh 
and  may  intersect  the  grids  in  an  arbitrary  fashion.  Our  method  is  based  on  an 
inverse  Lax-Wendroff  procedure  for  the  inflow  boundary  conditions.  Extensive  nu¬ 
merical  examples  are  provided  to  illustrate  that  our  method  is  high  order  accurate 
and  has  good  performance  when  applied  to  one  and  two  dimensional  scalar  or  system 
cases  with  the  physical  boundary  not  aligned  with  the  grids  and  with  various  bound¬ 
ary  conditions  including  the  solid  wall  boundary  condition.  In  [6],  this  technique  is 
generalized  to  compressible  inviscid  flows  involving  complex  moving  geometries.  Nu- 
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merical  examples  in  one  and  two  dimensions  show  that  our  boundary  treatment  is 
high  order  accurate  for  problems  with  smooth  solutions.  Our  method  also  performs 
well  for  problems  involving  interactions  between  shocks  and  moving  rigid  bodies.  In 
[7],  a  more  efficient  implementation  of  this  technique  is  introduced  and  the  scheme  is 
applied  to  very  demanding  test  cases  in  detonation  wave  propagation.  The  method¬ 
ology  in  this  line  of  work  is  expected  to  significantly  expand  the  applicability  of  high 
order  finite  difference  scheme  on  Cartesian  meshes  for  solving  problems  in  complex 
geometry  with  stable  and  accurate  boundary  performance. 

In  [2,3],  we  have  developed  several  types  of  conservative  high  order  semi-Lagrangian 
finite  difference  WENO  methods  for  advection  in  incompressible  flow  with  applica¬ 
tions  to  the  Vlasov  equation.  Each  type  has  its  own  advantages  and  limitations, 
however  collectively  these  methods  can  address  many  problems  in  plasma  physics 
very  efficiently.  In  [4],  we  have  developed  a  positivity-preserving  semi-Lagrangian 
discontinuous  Galerkin  with  applications  to  Vlasov-Poisson  systems. 

In  [9],  we  have  discussed  the  importance  of  well-balanced  schemes  for  moving- 
water  equilibria  of  the  shallow  water  equations. 

In  [10],  we  have  developed  a  WENO  scheme  with  subccll  resolution  for  computing 
nonconservative  Euler  equations  with  applications  to  one-dimensional  compressible 
two-medium  flows.  This  is  a  first  step  towards  the  development  of  a  robust  and 
accurate  solver  for  multi-medium  flows. 

In  [12],  we  have  systematically  studied  a  high  order  adaptive  finite  element  method 
for  solving  nonlinear  hyperbolic  conservation  laws. 

In  [13],  we  have  studied  a  strategy  to  improve  convergence  to  steady  state  solutions 
of  Euler  equations  with  high  order  WENO  finite  difference  schemes.  This  is  expected 
to  improve  significantly  the  efficiency  of  high  order  finite  difference  WENO  solvers 
for  steady  state  flows. 
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Journal  of  Computational  Physics,  v228  (2009),  pp. 8872-8891. 
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